Spin-Hall Conductivity in Electron-Phonon Coupled Systems 
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We derive the ac spin-Hall conductivity a s n(io) of two-dimensional spin-orbit coupled systems 
interacting with dispersionless phonons of frequency o>o- For the linear Rashba model we show that 
the electron-phonon contribution to the spin- vertex corrections breaks the universality of o~ s jj_(u>) 
at low- frequencies and provides a non-trivial renormalization of the interband resonance. On the 
contrary, in a generalized Rashba model for which the spin- vertex contributions are absent, the 
coupling to the phonons enters only through the self-energy, leaving the low frequency behavior of 
ct s h(w) unaffected by the electron-phonon interaction. 
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PACS numbers: 72.25.-b, 72.10.Di, 72.20.Dp 

The recent prediction of intrinsic spin currents gen- 
erated by applied electric fields in semiconductors with 
spin-orbit (SO) interaction Fj], nas attracted intensive 
research on the subject 0, y| encouraged also by poten- 
tial applications in spintronic-based devices. In such sys- 
tems, the spin-Hall conductivity ct s h = Jy /E x , where 
Jy z is a spin S z polarized current in the y-direction and 
E x is the electric field directed along x, arises from the 
SO dependent band structure which, for clean systems, 
leads for example to ct s h = — e/87r for a two-dimensional 
(2D) electron system with Rashba SO coupling .2] or to 
(t s h = — 3e/87r for a 2D hole semiconductor [(|. 

Of special interest for both applied and fundamental 
research is the role played by scattering events which 
have been shown to modify in an essential way the clean 
limit results. The most drastic effects are found in the 
2D linear Rashba model, where ct s h reduces to zero for 
arbitrarily weak impurity scattering 0, H, El E lilt - 
while the universal value <t s h(w) = —e/8ir is recovered 
for finite values of the ac field frequency ui in the range 
t -1 < ui < A 0,0], where r _1 is the impurity scatter- 
ing rate and A is the spin-orbit energy splitting. On the 
contrary, in 2D hole systems with weak (short ranged) 
impurity scatterirrg^Ts H (cj) remains equal to — 3e/87r for 
< uj < A 0, 0, E|> while it becomes dependent on 
the imp urity potential if this has long range character 
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So far, the study of scattering effects on the spin-Hall 
conductivity has been restricted to the case in which the 
source of scattering is the coupling of the charge carriers 
to some elastic impurity potential. This leaves aside the 
contributions from inelastic scattering such those pro- 
vided by the electron-phonon (el-ph) interaction which, 
in the materials of interest for the spin-Hall effect, ranges 
from the weak-coupling limit in GaAs [lfj to the strong- 
coupling regime in Bi(100) |l7| . 

Because of its dynamic and inelastic character, the el- 
ph interaction may affect the spin-Hall response in a way 



drastically different from static elastic impurity scatter- 
ing, questioning the general validity of the commonly ac- 
cepted forms of (7 s h(w) summarized above. Furthermore, 
the issue of the vertex corrections, which are responsi- 
ble for the vanishing of cr S H (lo = 0) in the impure 2D 
linear Rashba model 13 LUl EJ 113 ■ acquires a new 
importance, since these should be altered by the el-ph 
interaction. 

In this letter we report on our results on the spin-Hall 
conductivity <t s h(w) for 2D systems with SO interaction 
coupled with dispersionless phonons of frequency u>q. For 
a linear Rashba model we show that, in the frequency 
range r^ 1 < uj <C A (with A < ui ) where the universal 
value — e/87r has been predicted, the el-ph contribution to 
the vertex corrections reduces ct s h(w) to the nonuniversal 
value — e/[8-7r(l + A/2)], where A is the el-ph coupling 
constant. Furthermore, we find that the el-ph spin- vertex 
contributions renormalize also the interband transitions 
and provide a further reduction of <t s h(w) for lo > luo- 
On the contrary, in a 2D generalized Rashba model, for 
which the spin-vertex contributions are absent, the el-ph 
interaction provides only a trivial self-energy correction 
to the interband transition, leaving the low frequency 
part of <t s h(w) basically unaltered. 

We consider the el-ph interaction as given by the Hol- 
stein hamiltonian generalized to include SO coupling: 

k,a ka/3 

q qka 

where c\_ a and a q (ck Q and a q ) are the creation (an- 
nihilation) operators for an electron with momentum 
k = (k x ,k y ) and spin index a =t,l, and for a phonon 
with wavenumber q. et = h 2 k 2 /2m is the electron dis- 
persion, ujq is the phonon frequency and g is a momentum 
independent el-ph interaction (Holstein model) . The use 
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of a Holstein coupling permits to focus solely on the re- 
tardation and inelastic effects of phonons, disentangling 
the study from possible momentum dependences of the 
el-ph interaction. Furthermore, the Holstein coupling is 
partially justified, for example, by the results on surface 
states |l7j and by the reduced momentum dependence, 
compared to 3D electron gases, of 2D electrons coupled 
to bulk polar optical phonons |ljj. In the following, we 
shall also include the coupling to a short-ranged impu- 
rity potential V(r) = VS mp J2i — Ri)> where Ri are 
the random positions of the impurity scatterers. 

Let us start by considering a linear Rashba model, for 
which the SO vector potential is = jk(— sin</>, cos </>), 
where 7 is the SO coupling and </> is the polar angle. The 
electron Green's function of the interacting system is 

G(k, iuj n ) = \ ^ [1 + sf2 k • <r]G s (fc, iw„), (2) 

1 S=±l 

where f2k = (— cos tf>, sin </>) and G s (k, ioj n ) — [iu) n — El + 
H — £(iw n )] is the Green's function in the helicity ba- 
sis with dispersion Ef. — h 2 (k + sko) 2 /2m, ko = mj/h 2 
is the SO wavenumber, fi is the chemical potential and 
Lu n = (2n + l)nT is the fermionic Matsubara frequency 
at temperature T. Due to the momentum independence 
of g and Vj mp , the self-energy S(ia>„) is independent of k 
and reduces to 



E(iu3 n ) = 



W(iuj 



2Nn 



2 T, 

s=± 



kG s (k,iuj n i), 
(3) 



where iVo = m/lnfi 2 is the density of states per spin 
direction and 



W(iuj n - iui n >) = f"'" - A 

2-KTl 



{iu n - ki3 n t) - lo : 



2, (4) 



where r 1 = 2TrriiV i 2 np No is the impurity scattering rate 
and A = 2g 2 N$/(jjQ is the el-ph coupling. In writing 
Eas. l|3l4|) . we have employed the self-consistent Born ap- 
proximation for both impurity and el-ph scatterings. 

The equations defining the spin-Hall conductivity are 
obtained from the Kubo formula applied to the el-ph 
problem. Hence, the spin-current-charge-current corre- 
lation function is 

eh 2 ^f ^ — -V 

K{iv m ) = T2^r(iw i ,zw n )B 1 (it i ; i ,iw„), (5) 

n 

where v m — 2mirT is a bosonic Matsubara frequency, 
ui = uj n + v m and 

/dk 
—k 2, ^2sG- s (k,iuji)G s (k,iuj n ). (6) 
s 

The vertex function T appearing in Eq.JSJ) satisfies the 



following self-consistent equation 

TV- • x , . W ( iu} n> - iu n ) r ,. . . 

T(iuji,iuj n ) = 1 + ^2^ AN~k [B 2 (tujp,iu n i) 

n' ^ ^ 

+k B 3 (iu>i>,iu) n >)T(iu; ll ,iuj n ,)] , (7) 
where loi> = uo n i + v m and 

— k 2 ^sG s (k,iuji)G s (k,iLLi n ), (8) 

/dk 
— fc^ G a (k,iui)G a >(k,iu n ). (9) 

All integrations over the momenta k appearing in the 
above equations can be performed analytically, while the 
self-consistent equations (0 and J7J are solved numeri- 
cally by iteration in the Matsubara frequency space. Fi- 
nally, the (complex) spin-Hall conductivity 
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(10) 



is obtained from the retarded function K r {uj) = K(u> + 
iS) extracted from K(iv m ) [Eq.©] by applying the Pade 
method of numerical analytical continuation. Although 
our numerical calculations can be applied to arbitrary 
values of A/Ep, where Ep is the Fermi energy and A = 
2jkp is the SO splitting, the following discussion will 
be restricted to the weak SO coupling limit A/Ep <C 1, 
common to many materials, for which some analytical 
results can be obtained. In our calculations we have used 
T = O.Olwo (or T = O.OOlwo for the case shown in Fig. 
121, which is representative of the zero temperature case. 

We start our analysis by considering first the case 
uiq > A for which, as discussed below, the el-ph effects 
enter mainly through the real parts of the self-energy 
and of the vertex function. In Fig. ^ we show the real 
and imaginary parts of the spin-Hall conductivity for 
ujq = 1.5A and A = 0, 0.5, 1.0, and for weak impurity 
scattering 1/rA = 0.05. In the absence of el-ph interac- 
tion (A = 0), we recover the known results H,E3 charac- 
terized by the strong interband transitions at to — A and 
by the vanishing of a s n(uj) as lo — > 0. Furthermore, in the 
intermediate-frequency region 1/r < u <C A, Re<r S H(w) 
is almost ^-independent and matches the universal value 
— e/87r. This is better displayed in Fig. Efa) where 
the low-frequency behavior is plotted for 1/rA = 0.005. 
Upon enhancing A, two new features emerge. Namely, 
the frequency of the interband transitions get shifted at 
a lower (A dependent) value and, as also shown in Fig. 
Eta), the intermediate- frequency real spin-Hall conduc- 
tivity deviates from —e/8n, indicating that universality 
breaks down when A / 0. The origin of these features 
can be understood from the analysis of Eqs. (JSJ and 
0- In fact, at zero temperature and for A/Ep <C 1, the 
retarded function K R {u>) reduces to 0] 



eh 2 ^f 
4m 



r d€ rt 



+ +co,e-)B 1 (e++co,e-), (11) 
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FIG. 1: (color online) Real (a) and imaginary (b) parts of the 
spin-Hall conductivity for too > A obtained from the numer- 
ical analytical continuation of Eas.l(5jl- ll()|l (thick lines). The 
analytical formula (II 12 P is plotted with thin lines and is al- 
most indistinguishable from the numer ical resu lts. The peaks 
of Imtr S H(w) are centered at to/ A — \J\ + A/2/(l + A). 

where e± = e±iS. For lo < loq, the integration appearing 
in (|llfl restricts the e + lo and e variables to e + lo\ < loq 
and |e| < loq, for which the self-energy on the real axis 
can be well approximated by T,(x±) = — ARe(x±) =pi/2r, 
where X- — e_ and x+ = e+ + lo. In this way, the 
quite lengthy integral equation for T(e + + u>, e_), which 
can be derived from Eq.||7J) by following the method of 
analytic continuation described in Ref.|l8|], reduces to a 
simple w-dependent algebraic equation. Its solution for 
lo/A C 1 is T(cj) ~ lo/[(1 + X/2)lo + i/2r] and, since 
Bi(uj) is a constant for to /A -C 1, the low- frequency 
spin-Hall conductivity becomes a s n(uj) ~ — (e/87r)r(u>). 
We recover therefore the vanishing of a s n(u>) for lo 
while, contrary to the A = case, we find that 
Csh(w) is approximately equal to the nonuniversal con- 
stant -e/[87r(l + A/2)] for 1/r < lo < A. The break- 
down of universality at intermediate frequency reported 
in Figs, da) and Ufa) stems therefore from the el-ph 
contribution to the spin-vertex correction which, from 
Eas. l|7l9[) . governs the intraband contributions to <j s yi(lo). 
A more refined calculation which takes into account also 
the interband transitions leads to: 

e lo 

8-(l + |)-+^-[(l + A)^ + ^] 2 c 

(12) 

which is valid for lo < loq and arbitrary lo/A (for 
A/E F < 1). For A = 0, Eq.JTU is identical to the 
formula already published in Refs.@, H(ij . Instead, for 
A > we recover the intermediate frequency nonuni- 
versal behavior discussed above together with an el-ph 
renormalization effect to the interband transitions, which 
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FIG. 2: (color online) Low frequency behavior of a a u{w) for 
coo > A (a) and for u>o <C A (b). The thick lines are the 
numerical results, while the thin lines are Ea. l|12|l . In panel 
(a) they are barely distinguishable. 

now occur at a frequency lo — A* , where 

= yr+A/2 

1 + A y ' 

for 1/r A <C 1. When compared with the numerical re- 
sults of Figs, n and Et a ) ; equation H12fl is in excellent 
agreement for all frequencies lower than loq. As a mat- 
ter of fact, Ea. H12f> is in very good agreement with the 
numerical results also for lo > loq as long as loq > A 
while, for loq < A, the w-dependence of <7 s h(w) starts to 
be affected by the imaginary contributions of the el-ph 
self-energy and of the vertex function. These effects are 
visible in Fig. EJb), where we compare the numerical 
results for to — 0.05 A (thick lines) with Eq. ((T2l (thin 
lines). The deviation of Rea S H(w) from — e/[8rr(l + A/2)] 
for lo > loq stems from intraband transitions mediated by 
the phonons which, in analogy to the low temperature 
optical conductivity of the Holstein el-ph model 0, H3] , 
ensure conservation of energy and momentum. At higher 
frequencies, the real part of the el-ph self-energy goes to 
zero as Xloq/lo for large lo/loq, and the interband transi- 
tions occur at the unrenormalized frequency w«A. 

Having established that the nonuniversality of ct s h(w) 
at intermediate frequencies and the non-trivial renormal- 
ization (|13Jl have their origin in the el-ph contributions to 
the spin-vertex correction, we now turn to evaluate the 
el-ph effects when the spin- vertex corrections are absent. 
To investigate this point we have considered a 2D gener- 
alized Rashba model where the SO interaction is of the 
form fJ k = jk N (- sin N(p, cos N(/)) [H|]. For N = 1, the 
linear Rashba model discussed above is recovered, while 
for N = 3 this model describes a 2D hole gas subjected 
to an asymmetric confining potential. Because of the 
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FIG. 3: (color online) (a): Real part of the spin-Hall con- 
ductivity for a generalized Rashba model with N = 3. Thick 
lines are the numerical results while thin lines are Ea. 1)140 . 
(b): low-frequency region for ujq <C A. 



the 2D linear Rashba model we can make use of the 
equivalence between a s n(uj) and the longitudinal in-plane 
spin susceptibility HH> an d directly relate the 

poles of (Tsh(w) with the time evolution of the spin po- 
larization S y (t) |ld . which can be measured by vari- 
ous techniques |2lj. We find thus from Eq.JT^J that for 
rA 3> 1, S y (t) is a function oscillating with frequency 
A*, Ea. 1)13)1 . damped by an exponential decay with rate 
1/t s = 1/[2t(1 + A/2)]. On the contrary, the decay rate 
in the rA <C 1 limit is independent of the el-ph interac- 
tion at T = 0, and reduces to the Dyakonov-Perel value 

1/t s = A 2 t/2[H1. 
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angular dependence of fik for N ^ 1, the vertex correc- 
tions are absent Q, and the correlation function K{iv m ) 
is simply given by Eq.JSJ with T(iuJi,iu) n ) — 1 and with 
the prefactor multiplied by N . Furthermore, the func- 
tion Bifiui, iuj n ) is as given in Eq.JfjJ, with dkk 2 replaced 
by dkk 1+N and with dispersion El = h 2 k 2 /2m + sjk 3 . 
Contrary to the linear Rashba model, now all el-ph ef- 
fects arise solely from the self-energies contained in the 
interband bubble term B\. Hence, in the weak SO limit 
A/Ep -C 1, where now A = 2jkp, and by using the 
same approximation scheme as above, for A, uj < uiq the 
spin-Hall conductivity is easily found to be given by: 

eiV A 2 

CT sH (w) = — — -To 77Z J TT ,., 1a - (14) 

87T A^ — [(1 + A)cj + i/t\ z 

Contrary to Eq. l|12fl . the above expression predicts a low- 
frequency behavior unaffected by the el-ph interaction. 
Namely: ct s h(w) = —eN/8n for uj <C A. Furthermore, 
the interband transition frequency is renormalized only 
by the el-ph self-energy (mass enhancement) factor 1 + A: 
A* = A/(l + A), in contrast with Eq.JTSJ) where the el- 
ph contribution to the spin- vertex corrections contributes 
with a factor y/T+ A/2. This behavior is confirmed by 
our numerical results for N = 3 reported in Fig. |2la) 
(thick lines), which fully agree with Eq. 114)1 (thin lines). 
Furthermore, as shown in Fig. E2b) for ojq = 0.05A and 
1/rA = 0.005, for uj > ujq we find a weak deviation 
from Eq. 1)14)1 due solely to the imaginary part of the self- 
energy, in contrast to Fig. &h) where the spin-vertex 
corrections have a much stronger effect. 

Before concluding, it is worth discussing how our re- 
sults can be obtained experimentally. In particular for 
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